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Abstract 

We use gauge-gravity duality to model the crossover from a conformal critical point to a confin- 
ing Fermi liquid, driven by a change in fermion density. The short-distance conformal physics is 
represented by an anti-de Sitter geometry, which terminates into a confining state along the emer- 
gent spatial direction. The Luttinger relation, relating the area enclosed by the Fermi surfaces 
to the fermion density, is shown to follow from Gauss's Law for the bulk electric field. We argue 
that all low energy modes are consistent with Landau's Fermi liquid theory. An explicit solution 

is obtained for the Fermi liquid for the case of hard-wall boundary conditions in the infrared. 
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I. INTRODUCTION 



Much recent work |l|-l32l] has described the compressible, non-superfluid states of quantum 
matter by the methods of gauge-gravity duality. A variety of "non- Fermi liquid" and claimed 
"Fermi liquid" phases appear to have been obtained, but their correspondence with the 

ated phases of condensed matter models with short-range interactions 

are perhaps the closest to realizing 



known strongly-corre 
remains uncertain 



15 



29 



25j . The recent models of Refs. 
a conventional Fermi liquid state, in that they have N Fermi surfaces (with iV — > oo) 
enclosing areas which add up to the total fermion density, as required by the Luttinger 
relation. However, there are low energy excitations associated with infrared "Lifshitz" or 
related geometries, and these do not correspond to known low energy modes of a Fermi 
liquid. Moreover, there are also low energy excitations at small momenta associated with 
Fermi surfaces with vanishing Fermi momenta, and the extra low energy modes are possibly 
linked to each other. 

We will begin with conformally-invariant quantum critical point in 2+1 spacetime di- 
mensions (possible recent examples are in Ref. |33|, 134 ) . This has a familiar gravitational 
representation in a AdS4 geometry. Then we apply a chemical potential fi, and assume the 
system ultimately crosses over to a confining and gapless Fermi liquid state with a small 
number of Fermi surfaces (including the case with just a single Fermi surface). A Fermi 
liquid is a confining state because its only low energy excitations are Fermi surface quasi- 
particles which do not carry charges of the gauge-field of the boundary conformal field theory 



271 ]. Field-theoretic models which can exhibit such a crossover were discussed in Ref. 



27 



This crossover to confinement will be manifested in a deviation from the AdS4 geometry in 
the infrared. 

As an aside, we note that at non-zero density, it is possible for a state to be both confining 
and gapless even without the Goldstone bosons of a broken continuous symmetry, as in a 
Fermi liquid. On the other hand, at zero density, the usual assumptions of particle physicists 
apply, and the only possi ble g apless excitations in a confining state are Goldstone bosons. 

Initial works |l|, |2j, |4|t8|, lll| accounted for the chemical potential only via a Maxwell term 



for the bulk gauge field, and this led to a Reissner- Nordstrom black brane, with a AdS? x R 
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near-horizon infrared geometry at zero temperature. Subsequent works [9_|, 
included back-reaction on the metric from the matter in a Thomas-Fermi approximation, 
leading to a Lifshitz geometry in the infrared. However, it is clear that neither of these 
geometries can represent a true confining state (such as a Fermi liquid 271, |29[ ) , given the 
many emergent collective excitations in the infrared. 

At zero chemical potential, the crossover from conformality to confinement has previously 

into a "AdS 
for confining 



been understood via a termination of the anti-de Sitter geometry [35| 

soliton" geometry. Such a geometry has also been used more recently [37 

states at non-zero chemical potential in the presence of scalar fields. Here we will examine 

a terminated AdS4 geometry in the presence of a non-zero fermion density. For numerical 



simplicity, we will use a simple 'hard- wall' termination 39|, |40j, but we expect similar results 

It is important 



to apply to other confining boundary conditions in the infrared [41H44 
to note that while all excitations are gapped in a terminated geometry at zero chemical 
potential, this is no longer true once the chemical potential is taken large enough to obtain 
a non-zero density. Then there are gapless excitations across the Fermi surface. Indeed, 
these are the only gapless excitations, and this is the key to realizing a Fermi liquid. 

Another important feature of our analysis will be an exact treatment of the quantum 
mechanics of the fermions, avoiding the Thomas- Fermi limit. Our initial mean field theory 
in Section [Til will treat the bulk metric and gauge field classically, but we will discuss conse- 
quences of their quantum fluctuations in Section IIHI We describe cases where the boundary 
theory has a small number of Fermi surfaces, and focus on the case with a single Fermi 
surface. The previous Thomas-Fermi analyses effectively had a continuum of an infinite 
number of Fermi surfaces; equivalently, their fermion wavefunctions were delta-functions 
along the emergent spatial direction. The necessity of moving away from this limit has been 



previously noted [16, 



24. 281. 



The presence of quantum fermions implies that their contribution to the average local 
currents and stresses are not local functions of the background fields: they depend on the 
fermion eigenfunctions, which are controlled by the full spatial dependence of the back- 
ground fields. Consequently, the problem of determination of the background fields cannot 
be reduced to set of coupled ordinary differential equations. In this respect, we differ from 
other recent work 
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which approximates fermion bilinears with a zero frequency con- 
tribution. Our formalism focuses on determining the bulk eigenstates of the fermions, in the 
absence of external sources. Knowledge of these allows subsequent computation via spec- 
tral representations of the bulk-to-bulk Green's functions, and also of the boundary Green's 
functions, which represent the responses to sources. 



II. MEAN FIELD EQUATIONS AND THEIR SOLUTION 

We begin by setting up the framework for the simplest discussion of a Fermi liquid state. 

We will model a Fermi liquid in 2+1 dimensions by a theory of 'quantum electrodynamics' 
in a 4-dimensional spacetime with metric g. The emergent spatial-coordinate is z, and the 
metric is asymptotically AdS4 near the boundary z — > 0. We will work in Euclidean time, 
and assume the metric only has diagonal components non-zero. Apart from the terms which 
depend only upon g, the action is 



S 



d A x^fg 



^F ab F ab + % (^T M D M ^ + m^) 



(2.1) 



Here F ab is the gauge flux of a U(l) gauge field A^, which is the bulk representation of the 
globally conserved U(l) charge Q in the boundary theory. The bulk theory also has Dirac 



fermions with U(l) charge q and mass m > 0, and we will use the convention followed in 

n rn 

Refs. 16|, 123J for the Dirac matrices. 

In the present section we treat the gauge field A^ classically, while the fermions ip will be 
treated exactly. In the absence of sources, the only non-zero component of the gauge field 
is its time-component (the electrochemical potential) A t = i$(z), and we can write the free 
energy density as 



T 



^/wwg ->, 



F.D + m 



(2.2) 



where T is the temperature, and V is the spatial volume of the boundary theory. The 
determinant in Eq. (12. 2 p is to be evaluated as described in Ref. LzL Now our task is simply 
to solve the equations 

r -77 

, $(z ->•()) =/i, (2.3) 



for &(z), and describe the properties of the resulting state of matter. In particular, the 
charge density of the boundary theory, (Q) is given by 



(Q) 



dfx 



(2.4) 



Eqns. (12.31) and (12.41) imply that (Q) can be determined from Eq. ( 12.21) by Hamilton- Jacobi 
theory to be 

(Q) = -- 2 lim -f- \^gg* z g tt Hz)\ ; (2.5) 

e z z->o dz L J 

The right-hand-side is just the electric field along the z direction at the boundary, and this 

is Gauss's law for the "surface" charge density. 

We now show that the solution of Eq. (12 .3p is a manageable problem, and that the 
resulting solution obeys the Luttinger relation on the Fermi surface volume of a Fermi 
liquid. 

We make a simple choice of metric of a truncated AdS 4 space with unit radius 



ds< 



1 



— (dz 2 + dt 2 + dx 2 + dy 2 ) 



(2.6) 



which is terminated at a 



rard wall z = z m ; the same choice is made in a popular model of 
40| | . We will specify boundary conditions of the fields at z = z m 



confinement in QCD 39 
below. 

It is now useful to write down the explicit form of the Dirac equation, whose solution 
determines the determinant in Eq. (12. 2p . We follow Refs. |6|, |23| and reduce the 4-component 
Dirac equation to a 2-component equation; the equation with spatial wavevector k along the 



x direction and energy Eg(k) is obtained from Eqs. (3.11-3.13) in Ref. |23| to be 

( ia y — - a x — - ka z - q<S>(z) ) X e,k( z ) = E i( k ) Xe.k(z), 
\ dz z J 



(2.7) 



where x is a two- component spinor, a are Pauli matrices acting on the spinor space, and £ 
is a discrete label for the energy eigenvalues. These eigenvalues are quantized by our choice 
of boundary conditions. For z — > 0, we choose the usual solution in the absence of sources 



X(z -)• 0) ~ z n 



(2.8) 



For the boundary at z = z m , we require that the Dirac operator on the left-hand-side of 
Eq. ( 12. 71) is self-adjoint [45l-l48| . We normalize our eigenstates so that 



dz xlk( z )XeA z ) = !> 



(2.9) 



and then from the z derivative term in the Dirac operator, the self-adjoint condition is 

x\( z m)cr y X2(z m ) = 0, (2.10) 

for any two spinors Xi,2- We can achieve this by taking either their upper or lower compo- 
nents to vanish at z = z m , and we will choose the lower component in our numerical solution 
below. With the Dirac operator self-adjoint, the eigenvalues Ei(k) will all be real. 

We can now write down simple explicit expressions for the free energy and its functional 
derivatives. The fermionic contribution to Eq. (12.21) is obtained by filling up all the negative 
energy states at T = 0, and so the ground state energy density is 
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1 f Zm , fd^\ 2 v^ f d 2 k „„,„, „„ N , 



2e 2 



\dz J 



(2.11) 



The functional derivatives of E^{k) follow from the Hellman-Feynman theorem. Taking the 
functional derivative with respect to $(^), we have from Eq. (12. 3ft 



1 d 2 $ 
e 2 dz 2 



e 



d 2 k 

4T 2 



e(-E e (k)) X U(z)xe,k( z ) 



0. 



(2.12) 



This is nothing but Gauss's law for the bulk electric field, with the second-term equal to 
the bulk charge density: note that the latter is not simply a local function of $(z), unlike 



previous 



approaches [9j, 



13 



15 



19 



29j. We can now integrate Eq. (I2.12p over z from z = 



to z = z m , and obtain using Eqs. ( 12. 5 p and (12. 9p 

= ^E/S^(-^)) ( 2 - 13 ) 



1 d$ 



This equation is our Luttinger relation: the right-hand-side is the total area enclosed by the 
Fermi surface, and this equals (Q) modulo any electric flux leaking out from the IR boundary 



at z 



16 



25 



271 ]. as could happen in a non- Fermi liquid phase. We want a Fermi liquid 



solution here, and so we have our final boundary condition for the electrochemical potential 

= 0. (2.14) 



dz 



With this boundary condition, our solution obeys the conventional Luttinger theorem for a 
Fermi liquid. 

A. Numerical solution 

This solution will describe the solution of Eqns. ( 12. 7p and ( I2.12p . subject to the boundary 
conditions described above. 

First, at /i = 0, the Dirac equation is solved by Bessel functions 

m{z) oc V~z ( ~(k + E £ (k)) Jm+1 ^ MeZ) ) . (2.15) 

Using the boundary condition of a vanishing lower component at z = z m , the parameter 
Mi is determined by the £'th zero of the Bessel function of order m — 1/2, j m -i/2,£, with 
Me = j m -i/2,e/zm- The energy eigenvalues are 



E e {k) = ±^ k 2 + Ml (2.16) 

representing Dirac fermions of mass Mi in the confined state. All the negative energy states 
will be occupied in the Dirac sea, and we measure all charge densities and electrochemical 
potentials relative to this Dirac sea. The completeness relation for the wavefunctions Xi{k) 
implies that the charge density of the filled Dirac sea is z- independent; we explicitly verified 
the independence on z by a numerical sum over the normalized Eq. (I2.15P to a large value 
of £. This z-independence is associated with the z-independent measure in Eq. (I2.9p . and 
confirms our choice of boundary conditions as a consistent description of a gapped, confining 
state with the spectrum in Eq. (I2.16p . This (infinite) constant charge density is implicitly 
subtracted from all our expressions. 

Turning on a chemical potential, we see no change in the solution as long as \/a\ < Mi, 

6 



the smallest Dirac mass. So the Dirac vacuum is gapped and incompressible. As /j, crosses 
Mi, an additional band of previously positive energy states will be occupied, and we will 
obtain a gapless compressible Fermi liquid. To describe this state, we write 

Xe , k (z) = Cz m ( ^ ) (2.17) 

where C is a constant set by the normalization in Eq. (j2.9p . The wavefunction components 
obey the boundary conditions 

/i(0) = 1 
h{z m ) = 

Hz -+ 0) = ^— ^ z, (2.18) 

and the differential equations 



dz 

f: = _ {El (k) + q${z) - k)h{z) - —h{z). (2.19) 

dz z 

We made an initial (arbitrary) choice of $(2), and numerically solved the eigenvalue problem 
defined by Eqs. (I2.18P and (12.191) . After determining the eigenvalues and eigenf unctions, we 
computed the electric field, S(z), by integration 

±±0(-E £ (k)) dzxUz)xUz)- (2-20) 

71 Jz 

Determination of the electrochemical potential required one more integration 

$( z ) =fi- dzS{z) (2.21) 







It can be verified that the boundary conditions in Eqs. (12.31) and (12.141) are satisfied. After 
this determination of $(2), we returned to the solution of Eqs. (12.181) and (I2.19p . and iterated 
the procedure until the solution converged. We found little difficulty in the iteration, and 
the solution typically converged accurately after ~ 30 iterations. 

We show our results for the dispersion E^{k) in Fig. [TJ both with and without a chemical 
potential. The electric field and electrochemical potential computed from Eq. (I2.20|2.2ip 
are shown in Fig. [2} and the wavefunctions in Fig. |3j In computing the electric field, we only 
included the single 0. value associated with the shaded states in Fig. [H in the summation 
in Eq. (12.201) . Strictly speaking, we should also account for the fact that the states in the 



Ei(k) +qfi 




FIG. 1: Dispersion spectrum of the fermions. The blue lines correspond to the spectrum Eq. (|2.161 
at m = 1, z m = 3, and fi = 0. The red lines are at qfj, = 1.7 and q 2 e 2 = 3. The horizontal red line is 
at q/j,, and the shaded region shows the additional states filled by the chemical potential. Note that 
the two sets of band dispersions are offset even at large momenta and energies: this arises from a 
Hartree shift in the energies due to the added density of particles. The wavefunctions of the states 
at large momenta are not modified by this shift. At smaller momenta, both the k dependence and 
wavefunctions are different between the two sets. 



negative Dirac sea have also had their wavefunctions modified by the change in fermion 
density, and so will contribute to the electric field. However, we expect such contributions 
to be suppressed by the Dirac mass gap, and neglected them for simplicity. The more 
significant question is whether such corrections can be divergent after summing over the 
infinite number of occupied negative energy states, but this is difficult to answer reliably 
using numerics. From the perspective of the boundary theory, we expect that all ultraviolet 
divergencies are associated with the conformal field theory, and introducing the IR scale 
H does not introduce new UV divergencies; consequently we do not expect divergencies 
associated with the electric field. 

III. BEYOND MEAN-FIELD THEORY 

We now extend our results to include full quantum fluctuations of the quantum- 
electrodynamic theory in Eq. (12.11) . In principle, our arguments also allow for fluctuations 
of other fields, including the metric (suitably regularized). We will show that the state 
obtained in Section [III is a Landau Fermi liquid. 

The Gauss's Law result in Eq. (12. 12ft has a generalization in the full quantum theory. We 
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FIG. 2: Electric field and electrochemical potential for a non-zero fermion density. Parameters as 
in Fig. [TJ 



write the Ward identity associated with the equation of motion of $ as 

1 rf 2 ($) 



dz 2 



+ q — / —ImG R (u, k, z, z)n F (u) = 



(3.1) 



where G r (lj, k, z, z') is the exact bulk-to-bulk retarded Green's function for the fermion field, 
and np{ijj) = (e w / T + 1) _1 is the Fermi distribution function. In our mean-field theory, this 
Green's function has the spectral representation 



G R (u,k,z,z') = J2 



xj k ( z )XeA z ') 
u - E e (k) + W 



(3.2) 



and inserting this into Eq. ( I3.1i . we revert to Eq. ( 12.12p at T = 0. The expression in Eq. (13.2 




can be used to obtain response functions of the boundary theory by sending z, z' — y 

Now our central point is that we can treat the bulk G R as the Green's function of an 
inhomogenous 3-dimensional Fermi liquid with short-range interactions. This 3-dimensional 
system is confined in a potential along with z direction, and translationally invariant under 
the infinite x and y directions. The boundary conditions along the z direction ensure that 
the transverse components of the gauge field are gapped [40(| , and so there is no long-range 
interaction along the x and y directions. Thus, we can apply the conventional methods of 
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FIG. 3: Wavefunctions at k = Jzf with parameters as in Fig. Q3 The blue lines show the zero 
density state in Eq. (J2.15H . plotted at the wavevector with eigenenergy 1.7. The red lines are the 
state at the Fermi level in the finite density fermion system with \x = 1.7. The differences in the 
two lines are the consequences of the self-consistent electric field generated by the fermion density. 

many body theory to this inhomogenous system. These methods imply that the trace of the 
Green's function along the inhomogenous direction must obey a Luttinger theorem on the 
volumes of the Fermi surfaces associated with motion along the x and y directions. In other 
words, we can express the z motion in terms of a discrete set of "bands" with label £, and 
then the system looks very similar to a two-dimensional multiband solid state system. The 
Luttinger theorem applies to the sum over all bands, or equivalents, to the trace along the 
z co-ordinate. So we conclude that at T = 



dz 



du f d 2 k. 



7T 



4vr 2 



ImG H (u,k,z,z)6(-Lu) 



Areas enclosed by Fermi surfaces 

4^ ■ 



(3.3) 



The Luttinger result for the boundary theory now follows from application of Eq. ( 12. 5 p to 
the z integral of Eq. ( 13. ip . just as in the mean- field theory. 

Similar arguments, making the analogy to a multi-band two-dimensional Fermi liquid 
with short-range interactions, imply that this system realizes all other properties of a Lan- 
dau Fermi liquid. In particular, transverse gauge modes which are gapped and discrete 
at Gaussian order in the holographic theory [40(, will acquire Landau damping from the 
Fermi surface excitations at higher orders, and so their spectrum will eventually become 
continuous. Similar comments apply to the spectrum of gravitons and other bosonic modes. 
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IV. CONCLUSIONS 

We have presented a simple model for the holographic description of a Fermi liquid. We 



begin with a conformal field theory in 2+1 dimensions [33l . l34j . which has a holographic 
description by a theory on AdS4. Then we change its carrier density by applying a chemical 
potential \i. Using arguments that the Fermi liquid must be a confining state [273] , we 
assumed that the AdS 4 geometry was truncated by confinement at a distance z m along the 
holographic direction. The resulting state is described by 2 energy scales, /i and l/z m . We 
presented a theory for the crossover across these energy scales, and showed that a Fermi 
liquid is obtained at the lowest energy scales. The bulk gauge field was determined self- 
consistently, and this was closely linked to consistency with Luttinger's theorem. 

The main shortcoming of our analysis was that the confinement was essentially put in 
by hand, and not obtained by a self-consistent determination of the metric. In principle, 
confinement can be a direct consequence of introducing a non-zero \x on the conformal field 
theory. In this case, we expect the value of the confinement scale, l/z m , to determined by 
U, and to vanish as \i — > 0; such a situation would be analogous to 'deconfined criticality' 
[5ll . |52J. In particular, the scaling properties of the conformal field theory imply that z m = 
C/fi, where C is a universal constant characteristic of the conformal field theory. We leave 
the possible realization of such a scenario, and computation of C, as important problems for 
future research. 
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